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We test at the electroweak scale the recently proposed elaborate theoretical sce- 
nario for real-time dynamics of non-abelian gauge theories at high temperature. We 
see no sign of the predicted behavior. This indicates that perturbative concepts 
like color conductivity and Landau damping might be irrelevant at temperatures 
corresponding to the electroweak scale. 
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1 Introduction 



The study of SU(2) gauge theories at high temperatures is important for our un- 
derstanding of the electroweak theory in the early universe. One effect which has 
been of interest in this context is the baryon-number non-conservation in the elec- 
troweak theory, caused by the anomaly of the baryonic current. After the work 
of Kuzmin, Rubakov and Shaposhnikov |lj it was realized that the baryon num- 
ber non-conservation, caused by thermal field fluctuations between gauge-equivalent 
vacua with different winding numbers, could be large in the unbroken phase of the 
electroweak theory. A quantitative verification of the KRS-scenario requires non- 
perturbative real time simulations of hot thermal gauge theories, a task we still do 
not know how to perform from first principles. However, following suggestions by 
that one could use classical thermal gauge theory to address the question, it 
was shown in || that transitions between vacua with different winding numbers is 
unsuppressed at high temperatures in the unbroken phase of the electroweak theory. 
Starting with ||, ||, much work has since gone into refining the numerical tech- 
niques used and turning the qualitative statements into quantitative measurements 
II 0) S B OH 0) 01 an d also in understanding to what extent the notion of topology 



used in the continuum anomaly calculations is still valid in the lattice simulations 
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In order to address the determination of the rate of "sphaleron" transitions 
quantitatively one needs either to understand the corrections to the rate induced 
by using classical thermodynamics rather than the correct quantum field theory 
thermodynamics, or to develop better non-perturbative methods suited for real time 
simulations. The latter alternative is of course preferable since it might allow us to 
address many other non-perturbative questions which involve real-time dynamics of 
non-abelian gauge theories at high temperatures, by correctly incorporating thermal 
fluctuations in the ultraviolet. 

Much progress has been made in this direction over the last decade, starting 
with the concept of hard thermal loops, and culminating with the effective small- 
momentum, low-frequency theory of Bodeker ]TJ|] and its interpretation in terms 
of color conductivity jl5[ jiff] . In ordinary (abelian) plasma physics low frequency 
magnetic fields decay slower than naively expected due to Landau damping. It is 
now understood that the same mechanism applies to a pure non-abelian plasma, 
where the hard, high frequency modes couple to the low frequency magnetic modes 
much like the charged particles in an abelian plasma, and in this way produce a new 
time scale for the magnetic fluctuations. While the dominant long-range magnetic 
fluctuations will occur at the (non-perturbative) length scale of the order l/g 2 T, the 
lifetime of the fluctuations will be of the order u ~ g 4 T\n(l/g). 

On a more quantitative level the soft classical fields (momentum k < gT) couples 
to hard currents according to 

E = D x B - J hard , (1) 
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where 

J hard = &E + £, (2) 

and where the effective noise term £ is determined by the fluctuation-dissipation 
theorem: 

(ff (f , x)g(f , x')> = 2Ta 6vf*8(t - t')8(x - x'). (3) 

In (H) cr denotes the so-called color conductivity, which to a leading log approxima- 
tion is given by 

° = — , m 2 = \{9T)\ l = ^-9 2 THl/g). (4) 
7 3 lD7r 

m denotes the Debye screening mass and 7 the hard gauge fields damping rate (given 
here for pure SU(2) gauge theory). The derivation of ([I]) is valid in a frequency range 

a; <C A/y and k > g 2 T. (5) 

Strictly speaking, in order for (|5[) to be valid, one has formally to require that In 1/g 
be parametrically large. In this frequency range one can ignore the time derivative 
in (|1]). Then 

DxB = cxE + £. (6) 
Working in temporal gauge and ignoring non-linearities, one obtains from (^) 

k 2 A + ak = £, (7) 

which leads to the decay of gauge correlators: 

r k 2 -, 

(A(t)A(O)) ocexpj 1\. (8) 

The decay time r for k ~ g 2 T is thus of the order 

i.e., much longer than the non-perturbative length scale l/(g 2 T). 

While this picture no doubt gives an appropriate description of the long wave- 
length, low frequency thermal fluctuations for spatial scales up to 1/ (g 2 T) and time 
scales up to 1 / \g 4 T ln(l / g)) when ln(l/g) in some sense can be considered large, the 
question arises how well it describes "sphaleron" physics at the electroweak scale 
where a ~ 1/30, i.e. g ~ 0.65. 

Presently there are good indications that the picture works well at the above 
value of g. Real time computer simulations have been done, using various ways of 
implementing the hard currents || |8], |H]] The primary observable in the simulations 
has been the sphaleron rate. This is so for good reasons. Firstly, it is an impor- 
tant observable which might play a role in explaining the baryon asymmetry of the 
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universe. Secondly, it is believed to be dominated by long wavelength, classical 
thermal field fluctuations, i.e., one can hope that the classical theory, or at least the 
improvements, like Bodeker's effective field theory, will allow us to determine this 
(non-perturbative) rate by computer simulations. Until now all computer simula- 
tions basically agree, and there is a kind of consensus in the community that the 
sphaleron rate goes as predicted by Arnold, Son and Yaffe (ASY) and Bodeker. 

The purpose of this article is not to cast doubt on the effective theory of Bodeker 
and ASY, but only to point out that it is may not be clear that it can be applied to 
the electroweak theory for temperatures around the phase transition. The deriva- 
tion of the effective theory, in particular in the framework of ASY, uses perturbative 
concepts like real time gauge correlators and color conductivity, and the explicit be- 
havior of these objects as functions of momemtum k and frequency u are used. 
While some of these objects are gauge dependent they can nevertheless be mea- 
sured in the same computer simulations used to determine the sphaleron rate and 
one can directly check if their long wavelength, low frequency part show the be- 
havior predicted by the effective theory. If not, one could be tempted to conclude 
that one needs to go to higher temperature, and that one should not try to match 
the sphaleron rate to formulas based on the validity of the long wavelength, low 
frequency effective theory. 

The rest of this paper is organized as follows: in Section 2 we describe shortly our 
simulation setup and the set of gauge-covariant and gauge-invariant objects likely 
to be governed by a long wavelength effective theory. Our results are reported in 
Section 3. Section 4 concludes. 

2 Simulation setup and observables 

As mentioned in the Introduction, there presently exists no genuine non-perturbative 
way of simulating real-time processes in non-abelian gauge theories from first princi- 
ples. Solving the classical field equations of motion is as close as one can get to to a 
study from first principles, in the sense that no external parameters are introduced 
into the system. The caveat is, of course, that the classical thermal distribution 
is incorrect: while the thermal fluctuations in the quantum theory are cut off at 
scale T, the thermal fluctuation of the classical theory is only cut off at the cut-off 
scale A ~ 1/a, where a is the lattice spacing. The non-perturbative, long-distance 
magnetic sector of the thermal theory is not expected to be very sensitive to an 
exact location of the ultraviolet cutoff, as long as the cutoff is not too close to the 
magnetic scale. At the same time, the perturbative sector will be changed and, 
whenever one encounters the Debye mass one has to make a replacement: 

ml ~ g 2 T 2 -> g 2 T/a. (10) 

A typical high-temperature field configuration will be dominated by its large- 
momentum components, i.e., for the classical field theory the components of the 
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order of the cut off, for the full quantum theory the components of order T. If one 
is interested in the dynamic behavior of the system at large distances (from the 
magnetic screening length and above) one needs a way to filter out the prevalent 
large-momentum components of the field. One way to do so is by cooling. In 
the context of sphaleron transitions the method was first applied in the very first 



numerical work on the sphaleron transitions at the electroweak phase transition, [T? 
as well as in the study of the flow of eigenvalues of the Dirac operator in presence 
of lattice sphalerons [R|]. Later it was introduced as a tool to reduce the large 
momentum, short distance lattice artifacts of the sphalerons in real time simulations 
0. As discussed in ||, moderate cooling is best suited for this purpose, since it 
leads to an exponential decay of high-frequency modes. In this way the sphaleron 
profile, buried in the thermal fluctuations, will be enhanced. However, one can go 
a step further and simply view the cooling as a general procedure for integrating 
out the large momentum components of the thermal classical field theory, basically 
leaving us with the large distance physics of the classical theory. 

Our procedure is thus the following: first we generate a thermal field configura- 
tion for an SU(2) classical Yang-Mills theory on a lattice ||, for a given temperature 
T. We then let the system evolve according to the classical equations of motion. 
At selected instances along the classical trajectory we extract long-wavelength in- 
formation from the field configuration by cooling (we defer to the next section the 
discussion of how deep the cooling should be). Using the new configurations we can 
measure whatever observables we have in mind, and since we are probing the low 
momentum, classical sector of the theory, we should be able to compare our results 
with the prediction from effective theory given by ©-© provided we make the 
appropriate substitutions like ([TDD an d provided that we are in a coupling constant 
and temperature region where (|l|) and (§) are valid. 

The classical Hamiltonian dynamics is most conveniently studied in the temporal 
gauge where the electric field E(x, t) is the conjugate momentum to A(x, t). It is 
worth pointing out that this is also the gauge where eqs. ([![)- (|3[) naively make sense. 
As explained in |L6|, [L8| eqs. (j!])-® can be generalized to other gauges, but we will 



not need it here. We will study unequal-time correlators of the form 

C n (t) = ^Jd 3 x (O 1 (x,t)O 2 (x,0)), (11) 

where are designed to probe the long-wavelength properties of interest. In 
particular, if one is interested in the magnetic degrees of freedom, a natural choice 
would be the autocorrelators of the magnetic field tensor and of its covariant curl: 

1 (x,t) = 2 (x,t) = { D ^ ( ^ X ) ' t) (12) 

There is a special choice of 0\p that allows to determine the color conductivity a. 
Consider a correlator of B(x, 0) with (||) at time t. It follows that 

u ,_ Jd 3 x (DxB(x,0)-DxB(x,t)) 
a[) - Jd 3 x(DxB(x,0)-E(x,t)> [6) 
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if the time lag t is large compared to the autocorrelation time of the effective noise. 

On the other hand, we can also study the autocorrelator of the color charge 
density 

p = D E, (14) 

if we are interested in detecting soft longitudinal excitations of the system (plas- 
mons) . 

Note that all the space-local correlators of gauge-covariant quantities are invari- 
ant under the residual time independent gauge transformations which are not fixed 
by the choice of temporal gauge. However, if we transform these correlators away 
from the temporal gauge, they will in general depend on a Wilson line in the time di- 
rection, connecting x, to x, t. This being the case, two remarks are in order. First, 
it is unclear to us to what extent (|l3|) can serve as a genuine gauge invariant defini- 
tion of color conductivity, since it is based on correlations between gauge covariant 
objects. In fact, as emphasized in fl6| it is not very clear even in a perturbative 
framework how to define color conductivity beyond leading perturbative order. In 



16| it was defined to next to leading order simply as the coefficient which appears in 



the effective theory ([j])-([|), when formulated in temporal gauge or related, so called 
flow gauges, of which the Coulomb gauge is a limiting case. Our definition can be 
seen as an extrapolation of this philosophy: we have defined an effective theory by 
integrating out the high momentum part of the (classical) thermal theory and we 
then define a by ([[]) and @). Then a measurement of a(t) by (|i~3|) can be viewed as 
testing if our effective theory looks anything like (|1|) and (0). 

Secondly, it is not apriori clear what effect the connecting temporal Wilson line 
has on the characteristic time scale of correlators between gauge-covariant quanti- 
ties. For this reason, we also study the dynamics of a truly gauge-invariant object, 
5 2 (x, x). As our results show, there is hardly any difference in correlation time 
scales for gauge-covariant and for gauge-invariant quantities. 



3 Numerical results 

A real-time correlator C(t) of cooled classical gauge fields on a lattice depends on 
the following dimensional parameters: the temperature and the coupling constant 
in the unique combination g 2 T, the system size L, the cooling time r, and, finally, 
the lattice spacing a. As usual, we express all the other dimensional quantities in 
terms of the lattice spacing. The choice of other parameters is dictated by physical 
considerations. In particular, the inverse lattice temperature (3 = A/(g 2 Ta) is chosen 
within a range where the ratio of the perturbative Debye mass mo to g 2 T of the 
classical theory is close to that of the full SU(2) Yang-Mills theory at electroweak 
temperatures T ~ lOOGeV. Since we are interested in the dynamics of fields with 
momenta of the order of g 2 T, the dimensionless combination L/ (/3a) should be large 
enough in order to avoid finite-size effects. Most of our simulations were performed 
at L/(/3a) = 2.4. We verified that variations of L/((3a) around that value did not 
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have a measurable effect. Finally, the cooling time r should be large enough in order 
to suppress modes with momenta harder than those on the g 2 T scale. In most our 
simulations (g 2 T) 2 r = 3.84. We saw that further increase of r had little impact on 
the real-time behavior of the correlators. 

Our original motivation for measuring real-time correlators of cooled fields was 
to determine the color conductivity, as explained in the previous section. We will 
discuss our (thus far unsuccessful) attempt to determine a at the end of this section. 
However, the real-time correlators we measure are interesting in their own right, and 
we will first describe their properties as they transpired in the simulation. 

In Figures [I] through |^ we present the autocorrelators (D x B(x, t) ■ D x B(x, 0)), 
(B(x,t) -B(x,0)), and (£?(x, t)£?(x, 0)) - (5 2 (x, t)) (5 2 (x, 0)), respectively. Note 
that in the first two cases we determine autocorrelators of gauge-covariant quantities, 
whereas in the third case the quantity in question is gauge- invariant. In each case the 
autocorrelators are normalized by their value at the origin, while the time variable is 
expressed in units of 4/ {g 2 T). Notably, in all three cases the curves corresponding to 
different values of (3 coincide as long as the correlators retain a substantial fraction 
of their original value. This property is especially evident in the first two cases, 
where the error bars are smaller. 




Figure 1: The autocorrelator (D x B(x, t) ■ D x B(x, 0)) versus time t in units 
of (g 2 T) _1 for (3 = 8.33 (crosses), j3 = 10.0 (squares), (3 = 12.5 (triangles), and 
P = 15.0 (stars). 



To further quantify this behavior of the time scales, we introduce the integral 
autocorrelation time defined for an autocorrelator C(t) as 
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C(t)0 -6 - 
C(0) 




1.2 1.6 

ta/(3 
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Figure 2: The autocorrelator (B(x, t) ■ B(x, 0)) versus time t in units of (g 2 T) 1 for 
ft = 8.33 (crosses), /? = 10.0 (squares), ft = 12.5 (triangles), and ft = 15.0 (stars). 



where in our numerical estimates the upper limit of integration is replaced by a finite 
value t u , for which C(t u )/C(0) <C 1. In Figure f|we plot the dimensionless quantity 
4/(g 2 Ttj) as a function of 1/ft = g 2 Ta/4. Remarkably, in all three cases tj turns 

out to be of the order of g 2 T and shows little dependence on the lattice spacing 
throughout the range considered. There is therefore no evidence that in this range 
of the lattice spacings our cooled autocorrelators follow the ASY - Bodeker scenario, 
wherein the expected behavior is t j oc l/(g 4 T 2 a), up to logarithmic corrections. 

The autocorrelator (-Bf(x, t)B 2 (x, 0)) turned out to be a much noisier quantity 
than the other two, resulting in much larger error bars despite comparable sample 
sizes in all the three cases. Apart from this difference, the integral autocorrela- 
tion times behave very similarly for the gauge- covariant and for the gauge-invariant 
quantities, even though the gauge-invariant autocorrelator of the former necessarily 
involves a straight adjoint Wilson line connecting the points x, t and x, 0. 

Next, we consider the color charge autocorrelator (D • E(x, t)D ■ E(x, 0)). As 
shown in Figure |5|, this quantity is strikingly different from the magnetic-field au- 
tocorrelators considered earlier. The time scale for the color charge correlation is 
proportional to the lattice spacing and does not depend on g 2 T. 

This result can be contrasted with perturbative predictions. One would expect 
that the color-charge autocorrelator is dominated by the plasmon mode, whose fre- 
quency in the classical theory is of the order gJr /a and whose decay rate is of 
the order g 2 T. We observe none of these properties in the range of lattice spacings 
considered. 

Finally, we attempted to determine color conductivity a, as defined in the In- 
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C(t)0 -6 - 
C(0) 




Figure 3: The connected autocorrelator (Bf(x,t)Bf(x,0)} - (-B?(x, t))(B?(x, 0)) 
versus time i in units of (g^T) -1 for /3 = 8.33 (crosses), /3 = 10.0 (squares), /3 = 12.5 
(triangles), and (3 = 15.0 (stars). 



troduction, using Eq. [13]. As Figure || demonstrates, this attempt failed in two 
ways. First of all, a(t) does not appear to approach a constant for times in excess 
of the expected autocorrelation time of the noise term £ (and far in excess of the 
measured autocorrelation time of the noise). Secondly, the numerical value of a(t) 
is very small (less than 0.25/a) compared to the value expected in the ASY scenario 
(acr « 15). Given this small value of cr(t), it is not clear how neglecting the E term 
in Eq. [I] can be justified. 



4 Conclusion 

In summary, we see no sign of the ASY-B6deker scenario in the classical SU(2) 
theory, in the regime roughly corresponding to the electroweak scale. This negative 
finding does not rule out the ASY-B6deker scenario in general, but one is led to 
question its applicability outside the asymptotically weak-coupling regime, which in 
the classical lattice theory corresponds to asymptotically small lattice spacing. 

Neither our findings contradict earlier numerical data for the classical Yang-Mills 
theory, in particular, the sphaleron rate measurement by Moore and Rummukainen 
0. Their simulation was performed in the range of couplings (or lattice spacings) 
which overlaps the one considered here. Results of that work are consistent with the 
zero continuum limit of the rate, as predicted by ASY and by Bodeker. But they 
do not rule out a finite classical rate in the continuum. 

This brings us to the following methodological remark. The sphaleron transition 
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0.03 0.06 0.09 0.12 0.15 

1/(3 

Figure 4: Inverse integral autocorrelation time in units of g 2 T plotted against 
g 2 Ta for (D x B(x,t) • D x B(x,0)) (triangles), (B(x,t • B(x,0)) (crosses), and 
<5?(x,f)B?(x,0)> - (Bf(x,t)){B*(xM (squares). 

rate is a very important quantity in its own right, worthy of a careful numerical 
study. However, such study may not be the optimal way to test the theory of ASY 
and Bodeker. Precisely because the sphaleron rate is an essentially nonperturbative 
quantity, there is no easy way to disentangle the perturbative ASY-B6deker damping 
from genuinely nonperturbative effects. To compound the difficulty, topology on a 
lattice is ill-defined and requires special treatment. By contrast, a study like the 
one reported here attempts to make a direct contact with perturbation theory, by 
measuring objects such as color conductivity. It therefore may be better suited for 
testing perturbative predictions. 
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Figure 5: The autocorrelator (D • E(x, t)D ■ E(x, 0)) versus time t in lattice units 
for (3 = 10 (crosses), j3 = 12.0 (squares), and j3 = 14.0 (triangles) Note the complete 
or partial overlap of the data points. 
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